Quantum Nature of the Big Bang 
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Some long standing issues concerning the quantum nature of the big bang are resolved in the 
context of homogeneous isotropic models with a scalar field. Specifically, the known results on 
the resolution of the big bang singularity in loop quantum cosmology are significantly extended 
as follows: i) the scalar field is shown to serve as an internal clock, thereby providing a detailed 
realization of the 'emergent time' idea; ii) the physical Hilbert space, Dirac observables and semi- 
classical states are constructed rigorously; iii) the Hamiltonian constraint is solved numerically to 
show that the big bang is replaced by a big bounce. Thanks to the non-perturbative, background 
independent methods, unlike in other approaches the quantum evolution is deterministic across the 
deep Planck regime. 
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Thanks to the influx of observational data, recent years 
have witnessed enormous advances in our understanding 
of the early universe. To interpret the present data, it 
is sufficient to work in a regime in which space-time can 
be taken to be a smooth continuum as in general rela- 
tivity, setting aside fundamental questions involving the 
deep Planck regime. However, for a complete concep- 
tual understanding as well as interpretation of the future, 
more refined data, these long-standing issues will have to 
be faced squarely. Examples are: i) How close to the 
big-bang does the smooth space-time of general relativ- 
ity make sense? In particular, can one show from first 
principles that this approximation is valid at the onset 
of inflation? ii) Is the big-bang singularity naturally re- 
solved by quantum gravity? Or, is some external input 
such as a new principle or a boundary condition at the 
big bang essential? iii) Is the quantum evolution across 
the 'singularity' deterministic? In the Pre-big-bang and 
Ekpyrotic scenarios, for example, the answer has been in 
the negative Q . iv) If the singularity is resolved, what is 
on the 'other side'? Is there just a quantum foam far re- 
moved from any classical space-time, or, is there another 
large, classical universe? The purpose of this letter is 
to summarize results from recent analytical and numeri- 
cal investigations within loop quantum cosmology which 
address these and related issues. 

Loop quantum gravity (LQG) is a background inde- 
endent, non-perturbative approach to quantum gravity 
Loop quantum cosmology (LQC) focuses on sym- 
metry reduced models but carries out quantization by 
mimicking the constructions used in the full theory [3(. 
Results to date in this area fall in two broad categories: 
a) Resolution of the big-bang singularity using modifica- 
tions of the gravitational Hamiltonian due to quantum 
geometry (4|; and, b) Phenomenological predictions from 
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effective equations that incorporate the modifications of 
the matter Hamiltonians due to quantum geometry (see, 
e -g-i [MIS)' As in the first category, we focus on the more 
fundamental issues. While previous results showed that 
the LQC evolution does not break down at the singular- 
ity, as pointed out e.g. in 0, they did not shed light 
on what happened before. By constructing the missing 
conceptual and mathematical infra-structure, we show 
that the universe has a classical pre-big-bang branch, 
joined deterministically to the post-big-bang branch by 
the LQC evolution. Our detailed analysis of the Planck 
regime also provides tools to test the validity of assump- 
tions underlying phenomenological predictions. 

We will illustrate effects of quantum geometry on both 
the gravitational and matter Hamiltonians through a 
simple example: the spatially homogeneous, isotropic 
fc=0 cosmologies with a massless scalar field. Although 
the approach admits generalizations, we focus on these 
models because a singularity is unavoidable in their clas- 
sical theory. The question is if it is naturally tamed by 
quantum effects. The answer in the 'geometrodynam- 
ical' framework used in older cosmologies turns out to 
be in the negative |§]. For example, if one begins with 
a semi-classical state representing a classical universe at 
late times and evolves it back via the Wheeler DeWitt 
equation, one finds that it just follows the classical tra- 
jectory into the big bang singularity 0. In LQC, the 
situation is very different. This may seem surprising at 
first. For, the system has only a finite number of degrees 
of freedom and von Neumann's theorem assures us that, 
under appropriate assumptions, the resulting quantum 
mechanics is unique. However, for reasons we will now 
explain, LQC does turn out to be qualitatively different 
from the Wheeler-DeWitt theory |10| . 

Because of spatial homogeneity and isotropy, one can 
fix a fiducial (flat) triad °ef and its dual co-triad °oj l a . The 
SU(2) gravitational spin connection A l a used in LQG has 
only one component c which furthermore depends only on 
time; A l a = c °oj l a . Similarly, the momentum Ef canoni- 
cally conjugate to A\ — representing a (density weighted) 
triad — has a single component p; Ef ~ p (det °ui) °ef. 
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p is related to the scale factor a via a 2 = \p\. How- 
ever, p is not restricted to be positive; under p — > — p 
the metric remains unchanged but the spatial triad flips 
the orientation. The pair (c,p) is canonically conjugate: 
{c, p} = (87rG7/3), where 7 is the Barbcro-Immirzi pa- 
rameter. 

Quantization is carried out by closely mimicking the 
procedure used in full LQG jlOj . There, the elementary 
variables which have unambiguous operator analogs in 
quantum theory are the holonomies h of connections A l a 
and the (smeared) triads Ef. Now, background indepen- 
dence leads to a surprisingly strong result ^Jj : m essence, 
the basic operator algebra generated by holonomies and 
triads admits a unique irreducible, diffeomorphism co- 
variant representation. In this representation, there are 
operators h representing holonomies and E represent- 
ing (smeared) triads. But there are no operators rep- 
resenting connections A l a themselves. In the cosmologi- 
cal model now under consideration, holonomies along a 
straight line of (oriented) length fi (with respect to the 
fiducial triad °ef ) are almost periodic functions of c of the 
form N^c) := exp(ipc/2). (The are the analogs of 
spin- network functions in the full theory.) In the quan- 
tum theory, then, we are led to a representation in which 
operators N„ and p are well-defined, but there is no op- 
erator corresponding to the connection c itself (because 
the 1-parameter group iV M is not weakly continuous in 
fi). This new quantum mechanics is inequivalent to the 
Wheeler- De Witt theory already at a kinematical level. In 
particular, the gravitational part of the Hilbert space 
is now L 2 (^Bohr, rf/iBohr), the space of square integrable 
functions on the Bohr compactification of the real line, 
rather than the standard L 2 (R,d/j.) [1C|. While in the 
semi-classical regime LQC is well approximated by the 
Wheeler-DeWitt theory, important differences manifest 
themselves at the Planck scale. These are the hallmarks 
of quantum geometry 0, . 

The new representation also leads to a qualitative dif- 
ference in the structure of the Hamiltonian constraint 
operator: the gravitational part of the constraint is a 
difference operator rather than a differential operator as 
in the Wheeler-DeWitt theory. The derivation 0, [^J 
can be summarized briefly as follows. In the classical 
theory, the gravitational part of the constraint is given 
by / d 3 xe i ^> k e- 1 EfE^F abk where e = | detE] 1 / 2 and F* b 
is the curvature of the connection A l a . The par t of this 
operator involving triads can be quantized using a 
standard procedure introduced by Thiemann in the full 
theory. However, since there is no operator correspond- 
ing to the connection itself, one has to express F^ b as 
a limit of the holonomy around a loop divided by the 
area enclosed by the loop, as the area shrinks to zero. 
Now, quantum geometry tells us that the area operator 
has a minimum non-zero value, A, and in the quantum 
theory it is natural to shrink the loop only till it attains 
this minimum. There are two ways to implement this 
idea in detail. Here, we will use the one which has al- 
ready appeared in the literature 0, |jj [^j although it has 



certain drawbacks in the semi-classical regime, especially 
in more general models. The second method will be dis- 
cussed in the second of the detailed papers , which will 
also show that the quantum bounce and deterministic 
evolution across the Planck regime persist if the second 
and more satisfactory method is used. In both cases, it 
is the existence of the 'area gap' A that leads one to a 
difference equation. 

Let us represent states as functions 5'(/i, 0), where <f> is 
the scalar field and (modulo a fixed multiple of £ P1 ) the 
dimensionless real number [i is the eigenvalue of p |1C| . 
Then in LQC the (self-adjoint) Hamiltonian constraint is 
given by |!j 

= [S(M)]- 1 [C Pf (M)*(M + 4/* o ,^) + C°Oi)*(M,0) 
+ C~ (m)*(a» - 4/i , 0)] =: -e#( M , 4>) (1) 

where C+Qj) = \ |m + 3/x | 3/2 - \p + »o\ 3/2 I ; 

C-(p) = C+(fi - 4 Mo ); C°(fi) = -C+Qjl) - C"(/i) and 
where (6/8nj£ 2 Pl ) 3 / 2 B(n) are the eigenvalues of the oper- 
ator \p\~ 3 ' 2 j|. The fixed real number fj, is determined 
by the area gap; (8^7/6) /«o^pi = A. 

Now, in each classical solution, is a globally mono- 
tonic function of time and can therefore be taken as the 
dynamical variable representing an internal clock. In 
quantum theory, even on-shell, there is no space-time 
metric. But since the quantum constraint dictates 
how <f>) 'evolves' as <f> changes, it is convenient to 
regard the argument <f> in ^(/i, (j)) as 'emergent time' and 
fi as the physical degree of freedom. A complete set of 
Dirac observables is provided by the constant of motion 
Prf, and operators fj,\^ a determining the value of p, at the 
'instant' <fi = <f> . 

Physical states are the (suitably regular) solutions to 
Eq JT). The map n defined by n*(^» = *(-/!,</>) 
corresponds just to the flip of orientation of the spatial 
triad (under which geometry remains unchanged); n is 
thus a large gauge transformation on the space of solu- 
tions to Eq. (JTJ. One is therefore led to divide physical 
states into sectors, each providing an irreducible, unitary 
representation of this symmetry. As one would expect, 
physical considerations imply that we should consider the 
symmetric sector, with eigenvalue +1 of H [9(. 

To endow this space with the structure of a Hilbert 
space, we use the 'group averaging method' Qjj. The 
technical implementation of this procedure is greatly sim- 
plified by the fact that the difference operator O on 
the right side of (JIJ is independent of <f> and can be 
shown to be self-adjoint and positive definite (on the 
Hilbert space L 2 (-RBohr, B([i)dp,Bohr)- Since 8 is a dif- 
ference operator, the resulting physical Hilbert space H 
has sectors H e which are superselected; H = (B e T~te with 
e G [0, 2;U ]. States <f) in H e (are symmetric un- 
der the orientation inversion n and) have support on 
points fi = ±e + 4rt/i . Let us consider a generic H t . 
(The small technical differences in the exceptional cases 
are discussed in 0; they do not affect the main con- 
clusions.) Wave functions ^{fx,4>) solve ^ and are of 



positive frequency with respect to the 'internal time' <fi. 
Equivalently, they satisfy the 'positive frequency' square 
root of Eq (jTJ): 

-^0* = Ve* (2) 

and the inner-product is given by: 

(*l|^phy= E B0i)^M)9 2 M), (3) 

Ai£{±c+4/i Z} 

where as usual Z denotes the set of integers. On these 
states, the Dirac observables act in the expected fashion: 

p^ = -ihd^ W (4a) 
AL 0) - e^W-^) M o ) (4b) 

One can also begin with the complete set of Dirac observ- 
ables (0J and show that © is the unique inner product 
which makes them self-adjoint. 

To construct semi-classical states and for numerical 
simulations, it is convenient to express physical states 
as linear combinations of the eigenstates of p<p and O. 
We first note that, for \i ^> fj, a , there is a precise sense 
0, ^3 m which the difference operator 8 approaches the 
Wheeler DeWitt differential operator 9, given by 

(e/)(/i)-(16 7 rG/3) Ai 3 / 2 (^I/')'. (5) 

(Thus, if one ignores the quantum geometry effects, Eq 
reduces to the Wheeler-DeWitt equation d 2 ^ = 
— 0\I>.) The eigenfunctions 

e fc ( M ) = (1/47T) Itf/'e*^ (6) 

of are labelled by a real number k and its eigenval- 
ues are given by w 2 = (7rG/3)(16fc 2 + 1). The com- 
plete set of eigenfunctions efe(/x) of the discrete opera- 
tor O is also labelled by a real number k and e^(/i) are 
well-approximated by e fc (/i) for (i> fi 9] . The eigenval- 
ues cu 2 (k) of O increase monotonically with |fc|. Finally, 
the efc(/i) satisfy the standard orthonormality relations 
< ek\e' k >= 5(k, k'). A physical state <f>) can there- 
fore be expanded as: 

/oo 
dfctf(jfc) ei s) (/i) e"'* (7) 
-oo 

where ^(k) is arbitrary (but suitably regular), u(k) is 
positive and e^\n) = (l/\/2)(efc(/x) + e k (-[i)). Thus, 
each physical state is characterized by a free function 
^f(k). (For proofs and subtleties see [j|.) 

Since we have the explicit Hilbert space and a complete 
set of Dirac observables, we can now construct states 
which are semi-classical at late times — e.g., now — and 
evolve them numerically 'backward in time'. There are 
three natural constructions to implement this idea in de- 
tail, reflecting the freedom in the notion of semi-classical 
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FIG. 1: The absolute value of the wave function ^ is plotted 
as a function of <f> and jj, (whose values are shown in multiples 
of iio). For visualization clarity, only the values of greater 
than 10~ 4 are shown. Being a physical state, ^ is symmetric 
under fi — > — fj,. 



states. The main results in all three cases are the same 
. Here we will report on the results obtained using the 
strategy that brings out the contrast with the Wheeler 
DeWitt theory most sharply. 

As noted before, p^ is a constant of motion. For the 
semi-classical analysis, we are led to choose a large value 
p*^ (~S> h in the classical c = G = 1 units. In the closed 
model, for example, this condition is necessary to ensure 
that the universe can expand out to a macroscopic size.) 
Fix a point (/x*, (j) ) on the classical trajectory with p^ = 
p0 which starts out at the big bang and then expands, 
choosing fjb* 1. We want to construct a state which is 
peaked at (ft*, pi) at the initial 'time' (j> = 4> and follow 
its 'evolution' backward. Let \f(fc) be a Gaussian, peaked 
at a value k* given 

by P% = -(y/lGKGtP/3) k*. Set 

p CO 

*(m,^)=/ fi^laMe^ 1 '^-*') (8) 

J oo 

where <j>* = -^3/l6irG In |/x*| + (j> . It is easy to verify 
that 'J'fjU, <j) ) is the desired initial data, sharply peaked at 
P<t> = P% an d M = If evolved using the Wheeler-DeWitt 
analog —id^ — $ of Eq JSJ, it would remain sharply 
peaked at the chosen classical trajectory and simply fol- 
low it into the big-bang singularity |9j]. However, if it is 
evolved via @, the situation becomes qualitatively dif- 
ferent. The state remains sharply peaked at the classical 
trajectory till the matter density reaches a large critical 
value (which depends on pi), but then bounces, joining on 
to the 'past' portion of a trajectory which was classically 
headed towards the big crunch (see figures). 

To ensure that these results are robust, a variety of nu- 
merical simulations were performed to integrate Eq 
using the adaptive step, 4th order Runge-Kutta method. 
Due to space limitation, we will summarize only one of 
these. H ere, w e chose e = 2fi a and initial data with 
p*^ = 10 4 V Gh 2 , fi* = 10 5 /i o and the spread in the Gaus- 
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FIG. 2: The expectation values of Dirac observables p,\j, are 
plotted (in multiples of fj, ), together with their dispersions. 
They exhibit a quantum bounce which joins the contracting 
and expanding classical trajectories marked by fainter lines. 



sian &(k) given by Ap^/pJ = 7.5 x 10~ 3 , where Ap^ is 
the uncertainty in p^. The boundary of the numerical 
grid was chosen at 1.5/x* (where <j) )\ < I0~ 24 .) On 
the boundary Eq J5J was approximated by the Whceler- 
DeWitt equation and 'outgoing wave' boundary condi- 
tions were imposed. Results of the evolution exhibit a 
quantum bounce as shown in figures 1 and 2. Away from 
the Planck regime the uncertainties in the Dirac observ- 
ables are essentially constant. 
We conclude with a few remarks. 

1. The main limitation of this analysis is the restric- 
tion to homogeneity and isotropy. The approach can be 
readily extended to incorporate anisotropic models and 
potentials for scalar fields. However, incorporation of in- 
homogeneities has only just begun. The hope is that the 
deterministic evolution of LQC will enable one to evolve 
perturbations across the Planck regime. 

2. The dramatic difference between the predictions of 



the Wheeler DeWitt theory and LQC can be intuitively 
understood through effective equations which can be de- 
rived from Eq using certain approximations 9] . One 
finds that quantum geometry (which is ignored in the 
Wheeler DeWitt theory) modifies the Friedmann equa- 
tions. The modifications are significant only in the 
Planck regime and come with the sign required to make 
gravity repulsive. 

3. A common feature with the early LQC papers is that 
we did not have to introduce new physical input such as 
a boundary condition at the singularity. We only asked 
that the quantum state be semi-classical at late times. 
This is an observational fact rather than a new theoret- 
ical input or a philosophical preference. However, there 
are also notable differences from the existing LQC litera- 
ture. First, while much of the phenomenological work Q 
in LQC has incorporated quantum geometry effects only 
on the matter Hamiltonian, here they were incorporated 
also in the gravitational part. Second, we constructed 
the physical Hilbert space, Dirac observables and semi- 
classical states, thereby extracting physics of the Planck 
regime, going significantly beyond the demonstration of 
singularity resolution. Specifically, our results show that 
the quantum geometry in the Planck regime serves as a 
'quantum bridge' between large classical universes, one 
contracting and the other expanding. Finally, the idea 
that the scalar field can be used as an internal clock has 
appeared before, especially in However, that anal- 
ysis used conventional quantum mechanics rather than 
the Bohr compactification which descends from full LQC 
Our final physical Hilbert space is also different; it's con- 
struction is not motivated by the Kodama state. 
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